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Abstract 

Three local bifurcations in DC-DC converters are reviewed. They are period-doubling bifur- 
cation, saddle-node bifurcation, and Neimark bifurcation. A general sampled-data model is 
employed to study types of loss of stability of the nominal (periodic) solution and their connec- 
tion with local bifurcations. More accurate prediction of instability and bifurcation than using 
the averaging approach is obtained. Examples of bifurcations associated with instabilities in 
DC-DC converters are given. 

1 Introduction 

There have been many studies of instabilities of DC-DC converters [TJ [2j [31 [H [5] . Prom a practical 
perspective, it is useful to classify instabilities depending on how and in what range of operating 
conditions they arise. Bifurcation theory is a tool that facilitates the study of loss of stability 
and its implications for dynamical behavior. Upon loss of stability of a steady state solution of a 
dynamical system, typically a bifurcation occurs in which new steady states can arise. Thus, loss 
of stability of one steady state may lead to operation at a new steady state. A useful classification 
of bifurcations is that of local bifurcation vs. global bifurcation [6]. In a local bifurcation, the 
original steady state is an equilibrium point or limit cycle. In a global bifurcation, the original 
steady state has some other structure (say, an almost periodic solution, or a chaotic orbit). The 
paper focuses on local bifurcations. It is due to the fact, from a practical point of view, that 
these bifurcations can be expected to arise before any global bifurcation. Three typical local 
bifurcations of a periodic orbit are period-doubling (flip) bifurcation, saddle-node bifurcation, 
and Neimark bifurcation. The associated phenomena are subharmonic oscillation, jump, and 
quasi-periodicity, respectively. 

The most popular model for stability analysis of DC-DC converters has been the averaged 
model [7]. The averaged model generally makes the following approximations: 

1. The nominal steady state is an equilibrium. 

2. Equilibrium stability is of concern. 

3. Only one steady state is assumed. 



4. The duty cycle is a continuous-time variable and unbounded. 



For DC-DC converters, however, the realities are 

1. The nominal steady state is a periodic orbit, i.e., a limit cycle. 

2. Orbital stability is of concern. A periodic orbit is stable if any state trajectory starting 
close to the orbit moves towards the orbit. 

3. Many periodic and aperiodic steady states may coexist. 

4. The duty cycle is a discrete-time variable and bounded. 

Close to the onset of instability, the periodic nature of the steady state operating condition needs 
to be considered in order to obtain accurate results. Indeed, it has been reported [9] that 
averaging leads to erroneous conclusions regarding the onset of instability. Take subharmonic 
oscillations for example. Period-one and period-two orbits coexist and have the same averaged 
trajectories. Therefore, the averaged model can not distinguish the two orbits and determine 
their stability. 

This paper employs general sampled-data modeling [10], [TTJ [121 El E] of DC-DC converters. 
The bifurcations, indeed exists in real applications, are best explained by the sampled-data model 
instead of the averaged model. The sampled-data model does not assume those approximations 
mentioned above. Examples of instability in DC-DC converters are used to illustrate these 
bifurcations. 

The remainder of the paper is organized as follows. In Section [2J local bifurcations of a 
discrete-time system are summarized. In Section [3l a general model for DC-DC converters 
developed by the authors in [THrTS] is recalled. The three bifurcations are studied in Sections 0}- 
El Conclusions are collected in Section [71 

2 Local Bifurcations in Discrete-Time Systems 

In this section, the basic bifurcation theory used in the paper is recalled. For details, the reader 
is referred to [16) . 

Consider a discrete-time parameter-dependent system 

x n+ i = f(x n ,a), x n e~R N , qGR (1) 

The parameter a is called the bifurcation parameter. Suppose x = xq{q) is a fixed point of 
Eq. (P) for all a. Denote A(a) = f x (xo(a), a), the Jacobian of / with respect to x at (xo(a),a). 
The fixed point x = xq{o) is called a hyperbolic fixed point if A{a) has no eigenvalues on the 
unit circle in the complex plane. If a bifurcation occurs, then it must occur for a value a* of a 
for which A(a) is nonhyperbolic. There are three ways in which parameter variation can result 
in hyperbolicity being violated, and these are associated with three distinct bifurcations: 

1. Period-doubling bifurcation (the bifurcation associated with a real eigenvalue passing 
through the value —1): There is a curve of fixed point in the x-a plane on both sides of a = a* 
and a curve of period-two points on one side of a = a* intersecting with the first curve at 
a = a*. 

2. Saddle-node bifurcation (the bifurcation associated with a real eigenvalue reaching the 
value 1): There is a unique curve of fixed points in the x-a plane passing through (xo(a),a*) 
and locally lying on one side of a = a*. 

3. Neimark bifurcation (the bifurcation associated with a pair of complex conjugate eigen- 
values crossing the unit circle): There is a curve of fixed points in the x-a plane on both sides 
of a = a* and the emergence of a small-amplitude "invariant circle" around the fixed-point on 

one side of a = a*. 
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3 Sampled-Data Model for Voltage or Current Mode Control 



Without loss of generality, only continuous conduction mode is considered. A summary of the 
sampled-data modeling of closed-loop DC-DC converters discussed in p31 [15] is given. This 
model is applicable both to voltage mode control and current mode control. 

A block diagram model is shown in Fig. [TJ In the diagram, A\, € H NxN , B\,B<z € H Nx2 , 
C,Ei,E2 € R lxAr , and D € R lx2 are constant matrices, x G R , y £ R are the state and the 
feedback signal, respectively, and N is the state dimension, typically given by the number of 
energy storage elements in the converter. The source voltage is v s , and the output voltage is 
v Q . The notation v r denotes the reference signal, which could be a voltage or current reference. 
The signal h(t) is a T-periodic ramp with h(0) = Vi and h{T~) = Vh- In current mode control, 
it is used to model a slope-compensating ramp. The clock has the same switching frequency 
f s = 1/T as the ramp. Within a clock period, the dynamics is switched between the two stages 
S\ and S2- follows. The system is in S\ immediately following a clock pulse, and switches to S2 
at instants when y(t) = h(t). 

h(t) = V l + (V h -V l )(±modl) 
clock 

y = Cx + Du 



Switching 
Decision 



Switch to S*i or S2 



-1 


f ± 


= A\x + B\u 
= E x x 


S 2 : < 


X 


= A 2 x + B 2 u 




= E 2 x 



Figure 1: Block diagram model for DC-DC converter operation in continuous conduction mode 

Let x n = x(nT) and steady-state duty cycle D c = d/T. A periodic solution x°(t) in Fig. Q] 
corresponds to a fixed point x°(0) in the sampled-data dynamics. Using a hat " to denote small 
perturbations (e.g., x n = x n — x°(0)), The linearized sampled-data dynamics is 

x n+ i = <M n (2) 

where 

$ = e A ^ T - d Xl - (( A i- A ^°(d) + (Bi-B2)u)C Ald 
C(A lX °(d) + B lU ) - h{d) 

Local orbital stability of the converter is determined by the eigenvalues of denoted as a [$]. 
The periodic solution x°(t) is asymptotically orbitally stable if all of the eigenvalues of <!> are 
inside the unit circle of the complex plane. 

For comparison, the linearized dynamics of the state-space averaged model is 

i=(^ + (A '- A2)x r + f'- ft) " c )x ( 4) 
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where 



A ave := AonD c + Aoff(1 — D c ) 

Bzsre ■= -BonA: + -BoFf(1 ~ D c ) 
Aave — ^ave-^avef 

4 Period-Doubling Bifurcation: Subharmonic Oscillation and 
Eigenvalue Crossing -1 

In the period-doubling bifurcation, a 2T-periodic solution arises besides the original T-periodic 
solution. In most DC-DC converters, the period-doubling bifurcation is supercritical, where 
the 2T-periodic solution is stable and the original T-periodic solution becomes unstable. An 
illustration of such a bifurcation is shown in Fig. Throughout the paper, the stable solution 
is denoted as solid line; unstable solution, dashed line. 



Figure 2: Periodic solution before and after period-doubling bifurcation 




Consider the example [5] of a buck converter under voltage mode control shown in Fig. [3J 
Let T = 400 fj,s, L = 20 mH, C = 47 /*F, R = 22 VL, V r = 11.3 V, gi = 8.4, Vi = 3.8, V h = 8.2, 
(then h(t) = 3.8 + 4.4[y mod 1]), and let v s be the bifurcation parameter. 

Let the state be x = (iL,vc)' one has 
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The circuit undergoes a series of period-doubling bifurcations beginning at v s = 24.5. The 
eigenvalues of $ (i.e. 0"(<I ) )) as v s varies from 13.1 to 25.068 is shown in Fig.|H They are calculated 
from Eq. ([3|), while [T7] obtains the same graph by numerical estimation. One eigenvalue of $ 
is —1 when v s = 24.5, where the period-doubling bifurcation occurs. 

After period-doubling bifurcation, the original periodic solution becomes unstable, and a 
stable 2T-periodic solution arises. Take v s = 26, for example. The unstable T-periodic solution 
and the stable 2T-periodic solution are shown in Fig. [5j 

The averaged model of this example has been studied in [8] and found to be stable for v s 
varying from 15 to 40. Therefore the averaged model does not predict the period-doubling 
bifurcation accurately. 
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Figure 3: Buck converter under voltage mode control 




Figure 4: cr(<E>) as v s varies from 13.1 to 25.068 

5 Saddle-Node Bifurcation: Jump and Eigenvalue Crossing 1 

In the saddle-node bifurcation, a stable T-periodic solution collides with an unstable one at the 
bifurcation point, and no periodic solution exists after the bifurcation. This may explain some 
jump phenomena, and sudden appearance or disappearance of the nominal periodic solution in 
DC-DC converters. An illustration of such a bifurcation is shown in Fig. [6j 

Consider a buck converter with a discrete-time controller (Fig. [7]), where T = 400 (is, L=20 
mH, C = 47 (iF, R = 22 $7. Let v s be the bifurcation parameter and it is varied from 18.5 V to 
20.5 V. For duty cycle 0.7, the nominal inductor current is about I p = 0.6785 and the nominal 
output voltage is about V p = 14.0263. The switching decision in the cycle, t S [nT, (n + 1)T), is 
designed as follows (similar to a leading-edge modulation where the switch is off first and then 
on in a cycle) : the switch is turned off at t = nT and turned on at t = nT + d n . The switching 
instant d n is updated by d n = ^(0.3T — ki(i n — I p ) — k v (v n — V p )), where ki = —8.574 x 10~ 4 
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Figure 5: Unstable T-periodic solution (dashed line) and stable 2T-periodic solution for v s = 26 



and k v = 5.53 x 10 5 are feedback gains, and £ is a limiter (to limit the duty cycle within 1): 



The discrete-time control law above leads to new dynamics. It produces different static and 
periodic solutions for different v s . First, the switch being always on is a possible operation under 
some circumstances. When the switch is always on, d n = for any n. The steady state solutions 
are constant instead of being periodic: v {t) = v s and i{t) = v s /R. From Eq. ©, the following 
inequality needs to hold in order to make d n = 0: 

0.3T - ki(i n - I p ) - k v (v n - V p ) = 0.3T — — I p ) — k v (v s — V p ) 

k- 

= 0.3T + kilp + k v V p - (-^ + k v )v s 
< 

Therefore for v s > + k v )/(0.3T + kjp + k v V p ) = 19.213, the switch can be always on. 

However, the switch being always on is not the only possible operation for v s > 19.213. For 
v s € (19.213, 20), there are another two periodic solutions: one is stable, another one is unstable. 

Take v s = 19.9, for example. One stable periodic solution with duty cycle 0.6267 and one 
unstable periodic solution with duty cycle 0.7878 are shown as the solid line and dashed line 
respectively in Fig. [8l The stable one has output voltage around 12.5 V; unstable one has output 
voltage around 15.7 V. As v s is further increased, these two periodic solutions move closer and 
collide when v s = 20. For v s = 20, one eigenvalue of $ is 1 and a saddle-node bifurcation occurs. 
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Figure 6: Periodic solution before and after saddle-node bifurcation 
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Figure 7: Buck converter under discrete-time control 




Figure 8: Stable periodic orbit (solid line) and unstable periodic orbit (dashed line) for v s =19.9 



If v s is increased a little bit above 20, the operation suddenly jumps to the situation where the 
switch is always on and the output voltage jumps from 14 V to 20 V. 

The circuit is simulated for v s £ [18.5, 20.5] and the resulting bifurcation diagram is shown 
Fig. El In the figure, the upper solid line is for the stable static solution when the switch is 
always on (hence duty cycle is 1 and d n = 0); the dashed line and the lower solid line are for 
unstable and stable periodic solutions respectively with duty cycle less than 1. For v s below 
19.213, there is only one stable periodic solution and the output voltage is regulated below 11 
V. 

6 Neimark Bifurcation: Quasi-Periodicity and Eigenvalues Cross- 
ing Unit Circle 

When a Neimark bifurcation occurs, the original periodic solution (with frequency / s ) is mod- 
ulated by another frequency, ^ i.e., times the argument (i.e., phase) of the pair of 
eigenvalue of crossing the unit circle. The state trajectory will be periodic (phase-locking) if 
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Figure 9: Bifurcation diagram showing saddle-node bifurcation at v s = 20 



these two frequencies are commensurate; otherwise it will be quasiperiodic. 

Consider the example |18t p. 228] of a buck converter under voltage mode control shown 
in Fig. [TUl The system parameters are f s = 15 kHz, L=0.9 mH, C=22 /xF, R=20 v r =5 
V, Ri = #2=7.5 kfi, #3=60 Ml, C 2 =0.4 /iF, w s =30 V, V L =2.8 V, V v =8.2 V, (then = 
2.8 + 5.4[y mod 1]). All parasitic resistances are ignored. Here under the voltage mode control, 
the output voltage is sensed, scaled by 0.5, compared with v r = 5, and input into the error 
amplifier. Therefore, an output voltage around 10 V is expected. 
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The fixed point x°(0) is (0.25, 10, 0.39)'. The eigenvalues of $ are 0.8799 and 0.8797±0.4474i, 
which are inside the unit circle. Thus the periodic solution is asymptotically orbitally stable. 

As v s is increased from 30 V, the magnitude of the complex pair of eigenvalues begins to 
grow. For v s =36.9, the eigenvalues (0.8897 ± 0.4567i) exit the unit circle. Thus a Neimark 
bifurcation occurs. A low oscillating frequency ^Z(0.8897 + 0.4567i)= 1132 Hz modulating the 
original one f s is expected (Z denotes the angle in radian). Since these two frequencies are not 
commensurate, the steady state is quasiperiodic. 

For w s =30 (before bifurcation), the stable periodic solution x°(t) is shown as a solid line in 
Fig. [TTJ For -u s =50 (after bifurcation), the periodic solution becomes unstable (dashed line in 
Fig. Hip. A quasiperiodic state trajectory arises (Fig. 1 12 j) . coexisting with the unstable periodic 
solution. 

Output voltage waveforms of the quasiperiodic steady state and the unstable periodic so- 
lution are shown as solid line (with larger amplitude) and dashed line respectively in Fig. [T3] 
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for v s =50. From the figure, the quasiperiodic steady state has two oscillating frequencies as 
expected: f s modulated by a lower frequency around 1132 Hz. 

Next, the state-space averaged model (HJ) is analyzed. For v s = 36.9, a pair of eigenvalues, 
0.3 ± 7113. 5i, becomes unstable. The modulating frequency is 7113.5/27T = 1132 Hz, which also 
gives correct prediction. This is expected because the T-periodic solution can be averaged to 
an equilibrium and the quasiperiodic solution can be averaged to a low frequency (1132 Hz) 
oscillation. 




Figure 10: Buck converter under voltage mode control 
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Figure 11: Stable periodic orbit (solid line) for v s =30 becomes unstable (dashed line) for v s =50 



7 Conclusion 



Local bifurcations in DC-DC converters are studied using sampled-data models. The bifurcations 
considered are period-doubling bifurcation, saddle-node bifurcation, and Neimark bifurcation. 
Physical phenomena associated with these bifurcations are subharmonic oscillation, jump, and 
quasi-periodicity, respectively, and they are explained in details. Orbital stability is emphasized. 
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Figure 12: Quasiperiodic state trajectory in state space for v s =50 
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Figure 13: Quasiperiodic output (larger amplitude) and unstable T-periodic output for v s =50 
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Instabilities in DC-DC converters can be related to these bifurcations. These bifurcations, indeed 
exists in real applications, are well explained by the sampled-data model. The only bifurcation 
found predicted by the state-space averaged model is the Neimark bifurcation. 
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